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ABSTRACT 

Observation of the expansion history of the Universe allows exploration of the physical properties 
and energy density of the Universe's various constituents. Standardizable candles such as Type la 
supernovae remain one of the most promising and robust tools in this endeavor, by allowing for a 
direct measure of the luminosity distance-redshift curve, and thereby producing detailed studies of 
the dark energy responsible for the Universe's currently accelerating expansion. As such observations 
are pushed to higher redshifts, the observed flux is increasingly affected by gravitational lensing 
magnification due to intervening structure along the line-of-sight. We simulate and analyze the non- 
Gaussian probability distribution function of de/amplification due to lensing of standard candles, 
quantify the effect of a convolution over many independent sources (which acts to restore the intrinsic 
average (unlensed) luminosity due to flux conservation), and compute the additional uncertainty due 
to lensing on derived cosmological parameters. For example, the "degradation factor" due to lensing 
is a factor of three reduction in the effective number of usable supernovae at z = 1.5 (for sources 
with intrinsic flux dispersion of 10%). We also derive a useful expression for the effective increased 
dispersion in standard candles due to lensing, as a function of redshift. 

Subject headings: gravitational lensing — gravitation — cosmology: observations — cosmology: theory 
— supernovae 



1. INTRODUCTION 

The next generation of cosmological probes will explore 
properties of our Universe with unprecedented precision. 
The broad aim is to characterize the global nature of the 
cosmology through a mapping of the expansion history of 
the Universe, a{t). This history is fully encapsulated in 
the luminosity distance-redshift relation, which depends 
on global quantities such as the dimensionless dark mat- 
ter density, Q.m , dark energy density, Vlw , and parameters 
describing the precise nature of the dark energy (e.g. its 
equation of state). Measurements of the expansion his- 
tory thus enable determinations of these cosmological pa- 
rameters, and in particular shed light on the nature of 
the dark energy. 

Type la supernovae (SNe) are an essential tool in the 
direct mapping of the subtle changes in the distance- 
redshift relation due to the dark energy. Indeed, ob- 
servations of supernovae play a major role in our con- 
fidence that the current expan sion of t he Universe is 
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iTonrv et a l. 2003; Knop et al.ri200a: .Riess et all 200l] 
These (cali bratab le) sta ndard candles l|P hillips 1199,'^ 
iRiess. Press, fc Kirshnei<ll995t iWang et al.ll2003D give a 
particularly direct and simple translation from observ- 
able to expansion history: the redshift of a SN is a di- 
rect measure of the expansion factor a at emission, while 
the flux (or magnitude) of the SN relates to its distance 
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through a cosmological inverse square law, and hence 
to the lookback time t to the SN explosion. Thus an 
observation of a SN is very close to pure cosmography: 
drawing a map of the geometry of spacetime. 

We have strong indications from cosmic microwave 
background radiation measurements (in concert with 
large-scale stru cture measurements ) that the geometry 
of space is flat ijSpergel et al.ll200^ . corresponding to a 
critical density Universe with total dimensionless energy 
density = 1. This global statement breaks down on 
scales small compared with the cosmic (Hubble) scale, 
where a wealth of structure, from galaxies to clusters of 
galaxies, is apparent (and holds encoded within it a sep- 
arate brand of cosmological information). 

Large-scale structure also affects our goal of mapping 
the full spacetime geometry, representing the expansion 
history of the Universe. This structure, through its grav- 
itational potential, alters the propagation of light from 
distant sources, including standard candles, changing the 
received flux and hence compromising our measures of 
distance. Since structure on various scales fills the Uni- 
verse, every light ray is affected at some level by this 
gravitational lensing. However, since photons are con- 
served by this lensing, the mean flux over the sources is 
preserved. Gravitational lensing therefore adds a statis- 
tical noise to the cosmographic mapping, one that can be 
overcome by averaging over sufficient numbers of distant 
sources. 

In fJ21we discuss the properties of the lensing flux mag- 
nification distribution, and in ^Sjwe show how this affects 
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the distance-redshift relation. We carry this through 
in 21 to the influence on extraction of cosmological pa- 
rameters, yielding both a degradation in precision and a 
bias, and derive a prescription for the number of standard 
candles as a function of redshift required to provide av- 
eraging adequate to dilute the lensing degradation below 
a desired level. The conclusion in |S1 summarizes the cir- 
cumstances in which lensing has a significant impact on 
cosmological parameter determination. The Appendices 
discuss the non-Gaussian effects of translating between 
flux and magnitude measurements, and the influence of 
compact object lensing in a "clumpy Universe" picture. 

2. LENSING MAGNIFICATION DISTRIBUTIONS 

Consider a large number of perfect standard candles 
(sources with identical luminosities), distributed ran- 
domly on the sky at the same high redshift z. When 
observed in an inhomogeneous Universe, they will pos- 
sess a distribution of apparent luminosities, or mag- 
nitudes, due to alteration of the light ray bundle — 
gravitational lensing — by intervening structures along 
the various lines-of-sight. The form of these lensing am- 
plification probability distribution functions (PDFs) de- 
pends both upon the underlying cosmology (fim, 
dark energy equation of state ratio u'(z), etc.) and 
upon the nature of the structures causing the lensing 
(parametrized by cosmological quantities such as the 
mass clustering amplitude ctsi as well as details of the 
matter distributions: the mass function of galaxies, spe- 
cific density profiles of galaxies and clusters, etc.). 

To determine the lensing amplification PDF we uti- 
lize the Stochastic U niverse Method (SUM) presented in 
iHolz fc Waldl (CHi). This method is based upon a care- 
ful analysis of the assumptions needed for a physically 
reasonable model of a globally Robertson- Walker, but 
locally inhomogeneous. Universe. It calculates statistics 
for the full (weak and strong) lensing probability distri- 
bution of the observed image brightnesses of standard 
candle sources at any redshift, for a variety of cosmolog- 
ical models and a variety of mass distributions. It does 
this in a Monte Carlo fashion, randomly generating a 
portion of a Universe near a given photon line-of-sight, 
and calculating the lensing effects in this corner of the 
Universe. Repeating this many times, the full lensing 
distributions can be well approximated. 

Comparisons of SUM to conventional N-body ap - 
proaches can be found in lWang. Holz. fc Munshil l)2002|) . 
which also argues that the lensing probability distribu- 
tions have a universal form. This implies that fine de- 
tails of the distributions (e.g., those dependent upon un- 
known aspects of the matter profiles of the lenses) will 
have an inconsequential effect on broad characteristics 
of the PDFs, such as the variance, or the shift of the 
mode. As the current analysis is dominated by these 
broad features, fine details of the distributions will have 
little impact on our results. We use the SUM approach, 
ra ther than the universal (weak lensing) form presented 
in lWang. Holz. fc Munshil (|2002J, to ensure that any ef- 
fects due to the strong-lensing tail are properly mod- 
eled. Unless otherwise noted, all the results which fol- 
low are for a flat ACDM model, with f],„ = 0.28. The 
dark matter is taken to be smoothly distributed in ha- 
los (with density given by singular isothermal spheres or 
Navarro-Frenk- White profiles; details of the density pro- 



files d o not qualitatively change the results ijHolz fc Waldl 
119981) ). The further clumping of the dark matter into 
compact objects (e.g., MACHOs) will enhance t he ef- 
fects o f lensing (Amanullah 2003) (see also Goobar 'et alJ 
f200^ for the strong lensing case). 

The general features of a lensing probability distribu- 
tion function, P{fJ.), for observed flux /i, are as follows: 
P(/i) is peaked at a de-magnified value (compared to the 
"unlensed" pure Robertson- Walker result, which is nor- 
malized to fi = 1), with a long tail to high magnification. 
The mean of the distribution matches the Robertson- 
Walker value, (/i) = 1, due to conservation of photon 
number. As the mode is unequal to the mean, the dis- 
tributions are manifestly non-Gaussian. 

Most sources will be slightly demagnified by the pres- 
ence of structure, making them appear fainter and hence 
further away. If lensing were overlooked the resulting 
cosmography would be distorted, artificially enhancing 
the accelerating expansion caused by dark energy. This 
has a very small effect at low redshifts, however, since 
there is little optical depth to lensing (it is smaller by an 
order of magnitude than the detected acceleration l|Holj 
Il998|) 'l. At high redshifts it would continue to bias to- 
ward acceleration, diminishing the expected deceleration 
from matter domination. 

The faint demagnification of the majority of sources is 
compensated for by a small number of highly magnified 
sources. For sufficient statistics this high magnification 
tail is well sampled, and the mean of the distribution 
approaches the true source flux. Precision cosmology re- 
quires just such multiple measurements, automatically 
leading to the lensing distribution becoming better sam- 
pled and less biased, whether we explicitly recognize the 
presence of lensing or not. For high enough statistics, 
therefore, lensing magnification will not upset our in- 
terpretation of the acceleration and deceleration of the 
cosmic expansion, merely add noise to the precision mea- 
surements of its detailed behavior and the cosmological 
parameters derived therefrom. 

It is important to note that the lensing noise can 
be a signal in its own right. The statistical charac- 
teristics of the dispersion carry information about the 
matter power spectrum, e.g. ag (jFrieman. Holz. fc Kno^j 
120051) ■ as well as elucidatin g the nature of the dark 
matter lIMetcaff fc SilkllT999t ISeTiak fc HolzllT999l) . The 
more highly magnified sources can be useful cosmolog- 
ical probes as well, sensitive to both t he cosmography 
and the large scale structure prop e rties llKochaneld 20041: 

mm IHuterer fc Mari200l 
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such as supernovae allow greatly improved separation of 
the cosmological backgroun d informa tion and the lens 
model (Holz 2001: Qguri. Suto fc Turner..200a) . 

This paper concerns itself with the impact of lensing on 
determination of the distance-redshift relation, and the 
derived cosmological model. We quantify the meaning of 
sufficient statistics to moderate the two effects: the non- 
Gaussianity, in the form of the shift of the mode from 
the mean, causing a bias in the measured distance, and 
the dispersion of the lensing distribution, contributing 
an additional source of noise to the measurement, i i3.ll 
addresses the case of standard (or perfectly calibrated) 
candles, while i j3.2l treats the case of sources with some 
dispersion in their intrinsic luminosity. 2] then discusses 
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the implications for cosmological parameter determina- 
tion. 

3. SAFETY IN NUMBERS 

3.1. Perfect Standard Candles 

We begin with an ideahzed experiment: a large sample 
of perfect standard candles at high redshift. Without 
lensing, these would all be observed to have the same 
brightness: a delta function PDF (normalized at = 
1). We now add in the effects of gravitational lensing, 
which both contributes a width to the observed PDF, 
and shifts the mode of the PDF to slightly demagnified 
values. As already emphasized, the non-Gaussian lensing 
PDF preserves the mean: 



(1) 



This crucial property implies that, for sufficiently high 
numbers of observed high-redshift standard candles at a 
given z, the average brightness (in flux) will converge to 
the appropriate, unlensed brightness. It is to be noted, 
however, that the second moment of the lensing PDF 
doesn't necessarily converge. For the case of point-mass 
lenses, the probability at high magnification falls off as 
and so the contribution to the second moment is 
given by: 



(m^) = dfifi'^Pifi) (X / dfi/fi, 



(2) 



which diverges logarithmically at high magnification,^ 
emphasizing the non-Gaussian nature of the PDFs. 

It is to be expected that the effects of non-Gaussianity 
will be mitigated by observing sufficient numbers of SNe, 
and more fully sampling the lensing PDFs. With this in 
mind, we define PAr(/i) as the lensing magnification PDF 
for the mean magnification of a sample of N standard 
candles (at a fixed redshift). We calculate Pi(/i) via the 
SUM code. The distribution for higher numbers of stan- 
dard candles can then be calculated recursively: 



Pjv(/X)= dpdfi' PN-l{fi)Pl{fi')S 



iN-l)fl + 



N 



N dfjiPN-i{fji)Pi{N^ji-{N 



(4) 



This recursion becomes particularly straightforward us- 
ing spectral methods. Defining Pi(fc) as the Fourier 
transform of the lensing PDF Pi(/Lti), the convolution 
becomes 

Pn{p) = (27r)(^-2)/2iV / dfcPf (fc)e-*^-^^. (5) 



As expected, the convolution (eq. 01 or [SJ preserves the 
normalization and mean of the distribution, and the vari- 
ance shrinks as This formula reduces to a simple 
expression in the case of Gaussian or log-normal proba- 
bilities (see Appendix . 

^ In practice this is mitigated by effects such as finite source 
size and obscuration. If the second moment does diverge, then the 
distribution of observed brightnesses of large numbers of standard 
candles does not necessarily converge to a Gaussian distribution 
by the central limit theorem, and statistical intuition based upon 
normal statistics could lead us astray. 




Fig. 1. — Effective lensing magnification distributions for mul- 
tiple perfect standard candles, at 2: = 1.5 in a ACDM concor- 
dance cosmology. As more sources are observed the distribution 
approaches a Gaussian, and eventually converges on a (5-function 
at the unlensed magnification, = 1. 



The distributions at z = 1.5, for various values of N , 
are shown in Figure ^ Note that even for as many as 
50 SNe averaged together, the resulting distribution in 
magnification is still visibly asymmetric. This is a result 
of the high-magnification tail possessed by the lensing 
distributions. Figure |21 displays the shift of the mode of 
the distribution, as increasing numbers of SNe are ob- 
served. As expected, the curve asymptotes (slowly) to 
a value of 1, since the mean of the PDFs is always pre- 
served, and for large numbers of SNe the PDF should be 
well sampled. 

Figures 1^1 and 0] display the variance of the multiple SN 
lensing PDFs, as a function of the number of SNe. Since 
we are restricting our attention to smoothly clustered 
dark matter (as opposed to point masses like MACHOs; 
see [Amanullah (2003) for a treatment of that case), the 
variance for Pi(/i) remains well-defined and finite. The 
difference between a and the full width at half maxi- 
mum (FWHM) determination of the width is further ev- 
idence of the non-Gaussianity of the underlying lensing 
distribution: for a perfect Gaussian, a — FWHM/2.36. 
For = 1, the standard deviation is 1.61 times the 
FWHM/2.36. By A^ = 50, this factor has gone down 
to 1.35, indicating a more Gaussian-like distribution. 

3.2. Sources with intrinsic luminosity dispersion 

If the sources were perfectly calibrated candles, then 
the distribution of observed fiuxes would exactly mirror 
those shown in Figuren(i.e., refiect only the lensing mag- 
nification suffered during propagation). However, astro- 
nomical sources, even calibrated candles such as Type 
la supernovae, retain some intrinsic variation in their 
luminosity. This gives an innate a priori fuzziness in 
the distance-redshift relation, and the observed relation 
is thus a convolution of both the intrinsic and lensing 
fiux distributions. In other words, the distribution of 
observed fiux, F, is given by 

P{F)^ I dFo [ dMp'"*(Fo)p'"-(/x)(5(i^-FoM) (6) 



= / ^p'"*(^V'™^(A*), 



(7) 
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Fig. 2. — Mode (peak) of the lensing magnification distributions 
for multiple standard candles (see Fig. Q, as a function of the 
number of SNe averaged together. For a Gaussian, the mode equals 
the mean. The (very slow) convergence to the mean (at 1) for 
many sources averaged together can be seen. This is shown both 
for perfect standard candles ((Tint = 0), and for SNe with intrinsic 
noise given by a Gaussian in flux with (Tint = 0.1. This is for SNe 
at 2 = 1.5. 
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Fig. 3. — Width of the lensing magnification distributions for 
multiple standard candles (see Fig.0, as a function of the number 
of SNc averaged together. This is computed both as a variance, cr^, 
and as a full-width half-maximum, FWHM. For a Gaussian, one 
has (T = FWHM/2.36. The gap between the cr and FWHM curves 
is an indication of the non-Gaussianity of the lensing distributions. 
This is shown both for perfect standard candles (crjnt = 0), and for 
SNe with intrinsic noise given by a Gaussian in fiux with (Tint = 0.1. 
This is for SNe at z = 1.5. 



where Fq is the intrinsic source flux. In terms of mag- 
nitudes, the PDF for a single lensed, imperfect standard 
candle is given by 

pint+ions^^^ ^ y d^iP^"'%^i)P"''{m + 2.5logn). (8) 

To form the N-fold convolution of N lensed, imperfect 
standard candle SNe, one takes the lensing + intrinsic, 
single SN PDF of equation (0 as Fi , and plugs this into 
equation or (O. As before, the normalization and 
mean of Pi are preserved for Ppf. 

For the intrinsic supernova flux distribution we con- 
sider both a Gaussian in magnitude (which corresponds 
to a log normal in flux) and a Gaussian in flux. The SN 
distribution is characterized by a mean intrinsic magni- 
tude ((to(z)) = mo{z) or (niz)) = 1), and an intrinsic 
dispersion, crint(2^)- Canonically one uses crint(-z) = 0.1 or 
0.15 (constant with redshift). Figures show the re- 
sulting mode and variance, as a function of the number of 
observed SNe, for intrinsic noise given by a Gaussian in 
flux with standard deviation given by CTint =0.1. Note 
that the addition of intrinsic noise to the sources de- 
creases the shift in the mode of the magnification PDF, 
but increases the variance. 

3.3. Lensing as a function of redshift 

The previous section concerned itself with the distri- 
bution of the observed flux of multiple supernovae, av- 
eraged together. These distributions are of interest be- 
cause they emphasize that lensing preserves the mean 
(unlensed) flux, and thus for sufficient numbers of ob- 
served SNe the lensing can be averaged away. Figure^is 
a visual representation of the process of convergence to 
successively narrower Gaussians, approaching the delta- 
function limit. 

Guided by this convergence, we define an effective vari- 
ance, cr^ff, due to lensing. Consider some large number 
of standard candles, N, observed at a given redshift. For 
sufficiently high iV, the lensing distribution is found to 
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Fig. 4. — As in Fig.^l but on a log-log plot. The slopes of the 
standard deviation curves are very close to —1/2, indicating that 
the variance of the magnification distribution, including lensing, 
scales as l/N. 



be well-approximated by a Gaussian. If we designate 
the variance of this Gaussian as ct^ , we can then define 
the effective variance due to lensing, for a single standard 
candle, as cr^g = Naf^. Although the true distribution of 
a single standard candle is not well-described by a Gaus- 
sian of standard deviation a^ff, Figure ^ indicates that 
the scaling of the variance aj^ is roughly inverse linear in 
N. So (Toff gives an equivalent dispersion for each source 
when considering sufficiently large, but finite, numbers 
of standard candles. A useful rule of thumb for a reason- 
able use of (Toff is to require N greater than some lower 
limit calculated by demanding that ~95% of the prob- 
ability of the A'^-folded PDF be within 2a n this yields 
values of > 10. 

As SNe do not all occur at exactly the same redshift, 
we need to know over what redshift width we can apply 
this criterion. As we take SNe spread out in redshift, 
the lensing PDFs, and hence (Toff, vary. If we bound 
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this variation such that the 95% criterion above only di- 
lutes to a 95 ± 1% criterion, this is equivalent to using 
a 2 ± 0.1 tTAT allowance. So we may group SNe within a 
5% fractional deviation in cttv; since as we will see be- 
low, the standard deviation is linearly proportional to 
redshift, this imposes Az/z < 0.05. In other words, so 
long as there are N > 10 SNe in each redshift bin of 
width 0.1 z, the lensing can be effectively described by a 
Gaussian with standard deviation Ucff- (We emphasize 
that none of this necessitates any actual binning for the 
SN analysis.) For SNe that are more sparse in redshift, 
the full non-Gaussian lensing PDF should be used. Suf- 
ficient numbers allow lensing to be treated as a scaled 
noise added in quadrature; greater numbers of SNe are 
required to compensate for this noise — see SQl 

Lensing has little impact at low redshifts, where a short 
propagation distance implies a low optical depth. As the 
distance to the source increases, the lensing effects be- 
come more substantial. We thus expect aeS to increase 
with redshift. Figure plots as a function of z, for 
perfect standard candles, and for supernovae with intrin- 
sic noise given by a Gaussian in flux with CTint = 0.1. The 
results can be well approximated by a linear fit: 

CTcff' = 0.088Z, (9) 

for the case of perfect standard candles. In terms of 
magnitudes, we find 

aofF,m = 0.093Z (10) 
(to be used if taking (Tint as Gaussian in magnitude). 

4. COSMOLOGICAL PARAMETER ESTIMATION IN THE 
PRESENCE OF LENSING 

To map the expansion history of the Universe and re- 
veal properties of the dark energy, next generation ex- 
periments are being designed to observe the luminos- 
ity distance-redshift curve to one percent precision out 
to high redshift; for example, this is a primary sci- 
ence goal of the Supe rnova/ Acceleration Probe^ (SNAP; 
lAldering et al.l l2004f) . From the cosmological inverse 
square law the luminosity distance, d^, is related to the 
flux magnification, by /x ~ The errors on dis- 

tance are then given by dd-L/dL — (— l/2)5/i//i ~ 17^/2. 
It follows from ^(e.g., see Fig.EJ that the distance dis- 
persion due to gravitational lensing of any given SN at 
z > 0.5 is of order or greater than one percent. Figures^" 
01 illustrate the way to overcome this noise: observations 
of sufficient numbers of SNe per redshift "bin" convert 
a high-dispersion, non-Gaussian lensing distribution to 
a narrow, Gaussian-like one. For example, according to 
Figure 121 in order for lensing to contribute less than 1% 
distance uncertainty at z = 1.5 requires 70 SNe (with 
Cint — 0.1), as compared to 25 SNe if there were no lens- 
ing. 

It is to be noted that, since gravitational lensing is 
achromatic, the redshifts of the sources are unaffected. 
Thus the effects of lensing are confined to changing the 
inferred luminosity distance.'^ This compromises the 
determination of the luminosity distance-redshift curve, 
and thus impacts the estimation of cosmological param- 
eters. 

^ http://snap.lbl.gov 

^ Lensing amplification can also alter which sources enter a flux 
limited survey (Malmquist bias). However, most surveys are de- 
signed to be forgiving of small fluctuations in the flux threshold. 
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Fig. 5. — The effective standard deviation in flux (see text) is 
plotted as a function of redshift, including lensing. For perfect 
standard candles (dashed (red) line) this illustrates the pure lens- 
ing contribution; also plotted is the convolution with SN intrinsic 
dispersion given by a Gaussian with standard deviation 0.1 in flux 
(solid (black) line). 



4.1. Fisher approximation 

Given an observable quantity, and a model for its mea- 
surement errors, the most straightforward method for 
generating constraints on dependent cosmological pa- 
rameters is the Fisher matrix approach. For a super- 
nova survey, the data consists of magnitudes at several 
redshifts, and the Fisher matrix is 

_^ d{m){z)d{m){z) N{z) 

where N{z) is the number of supernovae in a redshift 
bin. This formalism can treat both offset effects in the 
observed magnitudes (rn(z)) and dispersion effects in am- 

We take as a set of cosmological parameters Oi — 
{Al, 57m, Wo, Wa}, where is a nuisance parameter in- 
volving the absolute magnitude of the supernova and the 
absolute distance scale (i.e. the Hubble constant), Q,m 
is the dimensionless matter density (we assume a flat 
Universe), and and Wa parametrize the dark energy 
equation of state ratio w{z) — wq + Wa{l — a). 

As a fiducial cosmology we adopt ri„j = 0.28, wo = —1, 
and Wa = 0, and either include a Gaussian prior of 
a{Vlm) = 0.03 or a WMAP cosmic microwave back- 
ground constraint on the distance to the last scattering 
surface. We also explore variations of the fiducial model. 

The increased magnification dispersion from lensing, 
as a function of redshift, is shown in Figure El and given 
in the fitting function of equation H10|) . The total disper- 
sion (J,n{z) is then obtained by adding in quadrature the 
intrinsic SN luminosity dispersion and the lensing effect: 

f^m(^) =f^cV™(2)+f^i'„f (12) 

One can verify from Figure [31 that the fitting form con- 
volved with the intrinsic dispersion is a good approxima- 
tion to the full, numerical solution. 

Before we study the effect of the increased dispersion 
on determination of cosmological parameters, we obtain 
a first indication of the influence of lensing by asking 
how many supernovae are required to offset the increased 
dispersion. We can then correct for this extra scatter 
by increasing the sample size. If we want the magnitude 
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Fig. 6. — Degradation factor due to gravitational lensing, as 
a function of redshift. This gives tlie fractional increase in the 
number of standard candles to be observed, to reduce the dispersion 
to what would have been observed without lensing. This is shown 
for two values of the SN intrinsic noise, given by Gaussians in flux. 



scatter to be at the same level as if there were no lensing, 
then we observe extra supernovae such that 



N{z) 



(13) 



From equations l(Tni) and l(T^ with CTint = 0.1, 0.15, 
we find that the excess numbers required are AN/Nq = 
{z/z^,y with = 1.08, 1.61 respectively (i.e., in the bin 
around half of the SNe are simply going to offset the 
increased lensing dispersion) . Figure illustrates this 
"survey bloat" , or conversely the degradation factor in 
the effectiveness of N SNe. At z = 1.5 one requires 
nearly three times as many SNe (with intrinsic dispersion 
fint = 0.1 in flux) to obtain the same dispersion as the 
case without lensing. In other words, each z = 1.5 SN 
has 1/3 the statistical weight it would have had in the 
absence of gravitational lensing. 

For sufficient number of SNe (given by the degrada- 
tion factor in eq. I13|l . however, lensing dispersion has 
no ill effect on cosmological parameter estimation. Fur- 
thermore, because the cosmological leverage occurs not 
merely from the highest redshift SNe, but from SNe at 
all intermediate redshifts as well, the relevant degrada- 
tion factor is much more modest than suggested in Fig- 
ure |S1 For a fiducial distribution of 40 low redshift SNc, 
10 per 0.1 bin at medium redshift, and 1 per 0.1 bin 
at high redshift (with no systematic uncertainties other 
than lensing), and adding a 0.03 prior on n„i, the effect of 
lensing from equation (|12|l is a 15-20% weakening of the 
constraints on wq and Wa- Furthermore, this can be un- 
done by increasing uniformly the number of SNe in each 
bin by a factor of 1.4, not the factor of three suggested 
by Figure for z = 1.5 alone. (These numbers hold 
with the further addition of a WMAP prior on distance 
to the last scattering surface. For a fiducial cosmology 
with SUGRA dynamical dark energy {wq = —0.82 and 
Wa = 0.58), rather than a cosmological constant, the 
weakening is 20-30%, though now a WMAP prior helps 
reduce it to 5-7%.) 

In addition to dispersion, lensing leads to a skewness 
in the convolved flux distribution. For a large sample 
this is not worrisome since lensing preserves the mean 
flux; a proper SN cosmology analysis treats the data in 
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Fig. 7. — The mode of the distributions, as a function of red- 
shift. This is shown both in magnitude (upper panel) and flux 
(lower panel), for either SN intrinsic noise (in magnitude or flux, 
respectively), o-jnt = 0.1 (solid black), or perfect standard candles 
(dashed red). 



flux (though it may be quoted in terms of magnitude) , 
and so there is no offset in the mean (see Appendix lAl 
for a discussion of bias induced by flux-magnitude con- 
fusion). But if the distribution is not well sampled then 
one might find a magnification offset from the mean. We 
characterize this by the mode of the distribution, which 
will be shifted to magnification lower than unity. Such 
an offset would affect cosmological parameter determina- 
tion, and we can calculate this within the Fisher matrix 
formalism. Figure |7| shows the mode offset, reasonably 
fit for a SN with CTint = 0.1 mag by 



Ato(z) = 0.025 (z- 0.7), z>0.7 



(14) 



and zero shift for z < 0.7. 

Such an offset propagates into a cosmological parame- 
ter bias via 

f dfTi 1 

{F-% dzN{z)Am{z) — ^y^, (15) 



AO, 



do. 



and gives a bias of less than 0.2ct (for the SUGRA case 
as well), where ct is the statistical error. Indeed, the shift 
of the mode is of even less concern in that equation H14|l 
refers to a single SN. For N SNe at a given redshift, 
we expect the magnification distribution to approach a 
Gaussian, and hence the mode to approach unity, with its 
deviation diminishing as I/a/ZV- This behavior is demon- 
strated in Figure [21 Since the Fisher matrix increases 
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proportionally with TV (it is also known as the informa- 
tion matrix), then from equation (|15|l we see that the 
parameter bias will scale as the mode offset, and thus 
vanishes as N~^^'^. For the number of SNe in current 
data sets, the lensing bias makes little difference (unlike 
the flux-magnitude asymmetry; see Appendix^) due to 
the large statistical error, but as the number of SNe in- 
creases the scaling will make the mode offset effect es- 
sentially negligible. 

The overall conclusion of the Fisher analysis is that 
lensing is a relatively benign error, ft is of significance 
only at high redshifts (z > 1), where there are cur- 
rently few SNe measurements, and so large statistical er- 
rors, and hence only modest cosmological leverage. Once 
many high-z SNe are observed, however, then lensing 
again is not a major source of error since the mean mag- 
nification averages to unity. Only when there are merely 
a handful of high redshift SNe per bin will lensing play 
a significant role. 

4.2. Monte- Carlo simulations 

While the Fisher analysis of the previous subsection 
gives a simple, intuitive feel for the effect of lensing on 
cosmological parameter determination, there remains the 
issue of the non-Gaussian nature of the gravitational 
lensing signal. A more robust, if computationally in- 
tensive, approach is to directly simulate data sets, and 
explore the ensuing scatter in best-fit parameter space. 
This presumes that the likelihood of a set of derived pa- 
rameters, given an observed data set, is equivalent to 
the probability of observing a data set, for a given set of 
parameters. 

We specify a cosmology, and calculate the lensing dis- 
tribution via the SUM approach. We then simulate an 
"observed" data set, with the SNe at pre-defined red- 
shifts, and the intrinsic and lensing noise for each SN 
drawn randomly from the appropriate distribution. For 
each data set we find the best-fit cosmology (using a max- 
imum likelihood approach). We repeat this procedure 
for many (^40,000) simulated data sets, and draw con- 
tours in parameter space encompassing the distribution 
of best-fit models (e.g. the la contour is the smallest 
contour which encompasses 68% of the simulated data 
sets). 

Figure IHl shows the resulting la and 2a contours, for 
a SN set roughly representing the current data (see the 
caption for details). The presence of lensing makes only 
a slight difference in the contours. This is because the 
data set under consideration has a small number of high- 
redshift (z > 1) SNe, and so the effects of lensing are 
mitigated. At 3cr, however, lensing has a more signifi- 
cant influence, greatly increasing the dispersion. In ad- 
dition, lensing due to compact objects (MACHOs, e.g. 
solar-mass black holes), will have a greater impact. We 
consider the extreme case of f/4 of the dark matter in 
compact objects, representing the maximum reasonable 
lensing effects for a given cosmology. The outermost 
(dotted) set of contours illustrate this extreme case; the 
additional lensing leads to a ~ 30% degradation of the 
parameter estimation (also see Appendix ^ on clumpy 
lensing). 

The moral to be drawn from this figure is that lensing 
is unlikely to have dramatic effects on the cosmological 
parameters estimated from current supernova data. A 




Wo 

Fig. 8. — Comparison of contours in parameter space, both with 
and without lensing. This is for observations of 20 SNe evenly dis- 
tributed in the range 0.01 < z < 0.1, 20 SNe in 0.1 < z < 0.5, 8 
SNe in 0.5 < z < 0.9, and 8 SNe in 1 < z < 1.7. The intrinsic 
dispersion of the supernovae is taken to be a Gaussian in magni- 
tude with standard deviation given by Uint =0.1 mag, with a prior 
on Vim given by Om = 0.28 ± 0.03. This data set is roughly com- 
parable to the current state of SN observations (which are statis- 
tically more numerous, but with greater individual errors). Solid 
(red) contours are without lensing. Dashed (green) contours in- 
clude lensing effects. Dotted (blue) contours represent an extreme 
lensing scenario, where 1/4 of the dark matter is in compact objects 
(MACHOs). 



similar conclusion holds if one adopts a CMB prior rather 
than a matter density prior. However, we expect increas- 
ing influence of lensing as observations attain higher red- 
shifts, but have not yet achieved sufflcient numbers of 
supernovae to compensate for the lensing degradation. 

Figure ini shows that once high redshift surveys do ob- 
serve many supernovae, the effective dispersion formula 
of equation Q provides good insight, and the Fisher 
matrix formalism, properly handled, can give a reason- 
able approximation. With 100 SNe at z = 1 and 100 at 
z — 1.5, the contour including lensing is bloated by 30% 
in Wo and 20% in Wa relative to a Universe without lens- 
ing. The Monte Carlo contour can be predicted by using 
the appropriate CToff, or equivalently degraded number of 
SNe. 

Specifically, we see that the lower right half of the 
Monte Carlo lensing contour (dashed (green) curves) 
matches the Monte Carlo contours generated without the 
presence of lensing, but with dispersion corresponding to 
(Teff added in quadrature (dotted (black) curves), or with 
merely intrinsic dispersion but numbers degraded by the 
appropriate factor from equation l|13() (dot-dashed (blue) 
curves) . The upper left half of the lensing contour is not 
as extended, and approaches a Monte Carlo contour gen- 
erated without lensing (solid (red) curves). 

We can understand this as follows: The models toward 
the lower right have equation of state ratios w < —1 
at redshifts z 1 (where lensing is important). For 
such values of u", the SNe should be dimmer than for our 
fiducial model. However, lensing magnification, /i > 1, 
can offset this dimming and make these models of nearly 
equal likelihood to the fiducial — i.e. they will lie within 
the confidence level contour. For models toward the up- 
per left, on the other hand, the opposite obtains: w > —1 
and ji <1. The /i > 1 part of the lensing PDF, convolved 
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Fig. 9. — At redshifts z > 1, lensing effects become more sig- 
nificant. Witli sufficient numbers of supernovae, however, the ef- 
fects can be treated via a Gaussian effective dispersion, or number 
degradation factor, and well approximated with the Fisher matrix 
formalism. Here contours represent parameter estimation from 40 
SNe at z < 0.1, 100 SNo at z = 1, and 100 at 2 = 1.5, includ- 
ing a prior on f2m given by Qrn, = 0.28 ± 0.03. The solid (red) 
contours are for SNe with intrinsic dispersion given by a Gaussian 
in magnitude with standard deviation (Tint = 0.1 mag, without 
any gravitational lensing effects considered. The dashed (green) 
contours are for the same data set, but with the effects of lens- 
ing included through Monte Carlo of th e PD F. The dotted (black) 
contour represents the CT^ff case (see where the same SN 

data set is considered, with the lensing effects approximated by an 
increased Ga ussian intrinsic dispersion in magnitude for the SNe 
(see eq. 1121 1. The dot -das hed (blue) contours demonstrate the 
degradation factor (see >I4.1I . This is for SNe with intrinsic noise 
given by Cint = 0.1 mag, but with the lensing approximate d by 
a diminution in the number of SNe, in accordance with eq. 1131 . 
yielding 54 SNe at z = 1 and 34 SNe at z = 1.5. As expected, the 
(Teff and "degradation" contours are essentially indistinguishable. 



with the intrinsic PDF, is well represented in its cosmo- 
logical effects by a Gaussian with a^ff. The fi < 1 part, 
even after convolution, is narrower due to the sharp low 
/X behavior of the lensing PDF, and so gives a Gaussian 
with dispersion closer to dint- This can be understood 
physically: gravitational lensing imposes a lower limit to 
the amount of rfemagnification, given by the empty-beam 
value (e.g., there is a minimum amount of focusing due 
to matter along any line-of-sight, given by the vacuum 
value). There is no similar limit to the magnification side. 
Thus, lensing is by its nature asymmetrical, and this is 
reflected in the contours. Note, however, that both sides 
can independently be reasonably approximated by the 
Fisher formalism with the appropriate dispersions. 

This general agreement is despite the fact that the PDF 
has not yet converged to a Gaussian form. Figure [HH 
clearly indicates that visible non-Gaussianity remains in 
the probability distribution even with 100 SNe at z = 
1.5: the lensing PDF is skewed from a Gaussian with 
the same variance (solid (black) curve vs. dashed (red) 
curve), and the variance is significantly different from 
the FWHM/2.36 (dotted (blue) curve is a Gaussian with 
the same FWHM as the lensing PDF). Thus, both the 
mode and the width of the magnification distribution 
exhibit non-Gaussianity — the central limit theorem still 
has a ways to go. Note, though, that the mean values for 
all the distributions agree (at the (unlensed) value of 1), 
and Figure O indicates that the cosmological parameter 
estimation is only dependent on broad characteristics of 
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Fig. 10. — When many high redshift sources are observed, the 
lensing magnification distribution will approach a Gaussian and 
eventually converge on a 5-function at the unlensed magnifica- 
tion, fj, = I. However, this convergence can be slow; the solid 
(black) curve shows the PDF for 100 SNe at z = 1.5. Resid- 
ual non-Gaussianities are clear: the lensing PDF shows deviations 
from both the dashed (red) curve representing a Gaussian with the 
same variance, and the dotted (blue) Gaussian distribution with 
the same full width at half maximum. However, we find that these 
residuals do not substantially affect cosmological parameter esti- 
mation, and it can be sufficient to account for lensing by means of 
a Gaussian effective dispersion cr^ff. 



the lensing PDF, sufficiently described by our acS- 

5. CONCLUSION 

Cosmological observations of the distance-redshift re- 
lation of Type la supernovae have revealed an extraor- 
dinary new component in the Universe: dark energy. To 
explore this in detail, we must comprehend the astro- 
physics that factors into the measurements and their in- 
terpretation. Realistic assessment of the data from large 
surveys over the next decade requires consideration of a 
realistic Universe. Structure on various scales fills the 
Universe, and every light ray is affected at some level by 
the resulting gravitational lensing. 

We have presented quantitative calculations of lens- 
ing effects on cosmological parameter estimation through 
observation of the luminosity distance-redshift relation. 
Using a Monte Carlo simulation code we derived the lens- 
ing magnification probability distribution function, and 
investigated its non-Gaussian characteristics, its behav- 
ior as the number of sources increased, and its redshift 
dependence. All of these have important ramifications 
for the interpretation of high-redshift observations. 

In particular, above redshift z ~ 1 lensing has a sig- 
nificant effect, degrading the statistical weight of sources 
by factors that exceed two (i.e. each source has less than 
1 /2 the cosmological leverage of one in a Universe with- 
out lensing), and increasing roughly as z"^. Standard can- 
dles at redshifts z > 2 will be severely impacted unless 
several hundreds can be observed at each redshift; this 
seems to doom hypothetical ultra-high-redshift standard 
candles such as gamma ray burst s or gravitational wave s 
from inspiraling black holes fsee lHughes fc Hol^ l)2003l) : 
iHolz fc Hughes! 1)2004 1) for a similar conclusion). 

We provide useful fitting formulas for the effective lu- 
minosity variance and mode offset induced by lensing, 
as a function of redshift, and compute their influences 
on cosmological parameter determination, checking these 
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against Monte Carlo simulations. Lensing effects, though 
still appreciably non-Gaussian even with 100 sources at 
a given redshift z > 1, can be well approximated by fit- 
ting functions for > 10 sources (within ^ O.lz). The 
fitting function for the flux dispersion induced by lensing 
is given by: 

cr^ff = 0.088Z. (16) 

From this expression, and from direct numerical calcula- 
tions, we find that TV « 70 sources are needed to achieve 
1% distance measurements at z — 1.5. This effective 
dispersion formalism, combined with the flux conserving 
properties of lensing (the mean magnification is unity). 



tells us that while lensing certainly cannot be ignored, it 
is tractable, with safety in numbers. 
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APPENDIX 

A. NON-GAUSSIANITY: LOG-NORMAL APPROXIMATION AND FLUX VS. MAGNITUDE 

While the full lensing PDFs derived from Monte Carlo simulations within the Stochastic Universe Method encapsulate 
the effects of structure on the luminosity distance-redshift relation, it is also useful to consider a simple approximation 
exhibiting non-Gaussianity. This will allow derivation of a number of analytic relations and help provide insight. 

Since astronomical observations of supernovae use a logarithmic flux system of magnitudes, with m = — 2.51ogj^Q F, 
it is natural to investigate a log normal distribution as a stand-in for the lensing magnification. This in fact has several 
attractive properties: a natural skewness, simple characterization, and Gaussian behavior near the peak. We do not 
claim it fits the derived lensing PDF, only that it has some related properties. 

A log normal probability distribution has the form 

p(^) ^ 1 J_e-("^''-^)'/(2S^), (Al) 



2^^S^l 

gC/+s /2^ with mean logarithm (In/x) = U, and variance cr^j — (/x^) — (/i)^ = e' 
that from the physical interpretation of the magnification, flux conservation ensures that (/i) = 1, and so for our case 
S'^/2. Thus we have a one parameter PDF with (In^) < 0, skewed negative as desired. 



and possesses mean (^) 



U 



1. Note 



We can now evaluate equation (jT)) and find that the convolution distribution for the observed flux is also a log normal 
distribution, broader and with different skewness: 



P(F) 



1 



1 



1 



-{\nF+r,io/b-Uf/[2{S^+af^Jb^)] 



(A2) 



where b — 2.5/ In 10 converts from flux in natural log to magnitudes, 
logarithmic flux is 

{\nF)^-^-^mo + U^{lnFo)- 



The mean flux is Fq, and the mean of the 



For observations using magnitudes, it is convenient to convert to the distribution in magnitude, m: 

1 



P{m) 



'I'K (7,1 



-(m-mo+f)(7)V(2o-„) 



(A3) 



(A4) 



{m){z) = mo(z) — bU{z) 
al{z) = al,-2b'U{z). 



where cr^{z) = crf^^ + 6^5*^. Also, recall that from flux conservation U ~ — S'^/2. Thus these probability distributions 
for source properties and lensing effects lead to a particularly simple result for the observed magnitude-redshift relation: 
there is an increased dispersion Um, and an offset in the mean. Explicitly, 

(A5) 
(A6) 

Recall that /7 < 0. The increased dispersion will weaken cosmological parameter determination, and the offset in the 
mean will bias the derived parameters. The magnitude of the effect is controlled by the quantity U = (In^). For 
example, the variance is doubled, relative to a 0.1 mag intrinsic dispersion, for U — 0.004. 

It is worth noting that no actual lensing is required to generate bias in the cosmological parameter determination. 
Image intensity can be expressed in either flux or magnitude, and the transformation between these systems of mea- 
surement gives a similar form of non-Gaussianity — that of the log normal distribution — with the attendant dispersion 
and bias. By confusing flux and magnitude values in an analysis of SN data, it is possible to severely com promise the 
values of the resulting cosmological parameters. Observing sources given by the numbers and redshifts in lRiess et alJ 
l)2004|) . with an intrinsically Gaussian flux distribution, but interpreting them in terms of magnitudes, turns a true flat 
cosmology with 0,^ — 0.28 and a cosmological constant into an apparently substantially closed cosmology [0,^ ~ 0.2). 
This might also be interpreted as a flat model with a dark energy different from wq = — 1, = 0, in particular 
wq < —1, Wa > . To avoid bias one must therefore carry out both the observations and the analysis in a fully 
consistent manner iHolz. Livio. fc Riesall200l . 
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B. CLUMPY UNIVERSE MODEL 

The central problem of incorporating lensing effects from structure in the Universe into the distance-redshift relation 
is calculation of the magnification probability distribution function (see |2]of the main text). The simplest distribution 
of lensing magnificatio ns, however, is no distribution at all. This is the approach of the clumpy Universe model of 
IDver fc Roedeil l)1973j) . This model does not ignore lensing, but rather assumes a special set of lines of sight — those 
that avoid mass concentrations. This preferentially avoids magnifications greater than unity. To conserve mass the 
model partially evacuates most of space, leaving an average density ap, where a is the smoothness parameter and p the 
globally averaged energy density, and concentrates the remaining mass into compact clumps that the "chosen" light 
rays never encounter. Furthermore, shear effects from these inhomogeneities are ignored. Thus the only light paths 
the observer deals with lead to demagnification (/i < 1) of the source, and overall fiux is not conserved (we mention a 
way around this below). 

We characterize this as shunning lensing rather than ignoring it, willfully neglecting data from sources touched 
by magnification; in general this is an even worse idea than completely ignoring lensing, due to the bias and flux 
nonconservation. However, we can use it as a simple, extreme case to examine the effect of lensing on the distance- 
redshift relation, and the ensuing cosmological parameter determination and bias (for example, see H4.2I and Fig. IHJ. 
Thus it is worth investigating some aspects of it here. 

The distance-redshift relation is no longer an integral over the Hubble parameter, but instead becomes a second 
order d i fferential equation f ollowing fro m the o ptical scalar equations. This was given for general equations of state by 
iKavseij ijunnublished. 1985^ and Lind'e^ l)1988D . and can treat inhomogeneities not only in the matter but hypothetical 
dumpiness in the dark energy component. The effect of dumpiness enters the Ricci focusing as [1 -|- w(z)][l — 
a.ujiz)]^.uj{z), with a^i giving the ratio of smoothly distributed to total energy density in a component characterized 
by equation of state w. 

We note several points: (1) a pure cosmological constant has no such Ricci focusing, (2) the effect is diminished as 
w — > — 1 but is enhanced for larger Q^, and (3) an evolving dumpiness can be treated straightforwardly (as detailed 
in iLinderl l)1988t) ). One could certainly assume the canonical smoothness a^, = 1 for a component other than matter. 
However we can also explore more arcane models, such as postulating that dark energy couples to (dark) matter, or is 
affected through backreaction on the expansion from nonlinear structure formation in the matter. Thus we examine 
both ttu) = 1 and = arm where am is the dumpiness of the matter component. Since the clumpy Universe model 
generally counts larger scale structure such as galaxies and clusters as a smooth component, only excluding compact 
objects with strong magnifications, we take the conservative position that all dark matter is smooth, and about half 
of the baryons as well, so am = 0.9. The effects will be greater for more extreme cases. If we take as a dark energy 
model the SUGRA model with wq — —0.82 and Wa = 0.58, then at z = the Ricci focusing due to clumpy dark 
energy would be roughly half that due to clumpy matter. Continuing to be conservative, we henceforth ignore that 
contribution and assume the dark energy is smooth. 

Before investigating the influe nce of clumpy lensing on cosmolo gical parameter determination, let us consider pa- 
rameter degeneracies. Recently, iCaldwell fc KamionkowskH l)2004l) pointed out that formally one could measure the 
curvature of space through the distance-redshift relation, with the curvature entering at third order in a low redshift 
expansion; however, the space geometry is entangled with the spacetime geometry, i.e. the expansion history. This 
usefully highlights the same point made in .Linder (1988) . Here we extend that result to show a six- fold degeneracy at 
0{z^) of the angular diameter distance: 

r{z) w z - (l/2)(3 + qo)z^ + z^[2 + (S1t/4)(2 + 5w - 3^ - - a - wa) + {l/2)ql] + 0{z^). (Bl) 

The notation X — J2w ■^^■wW/^t, where fix is the total dimensionless energy density, related to the curvature 
density by Ofc = 1 — f^T and to the deceleration parameter by go = f^T(l + 37Z))/2. Thus one must know the average of 
the equation of state ratio w, of the time variation w' = dw/dz, smoothness a, etc. to determine the spatial curvature. 
This points up the well known necessity for distance observations at many redshifts over a long baseline, e.g. z w — 2, 
to break the degeneracies and determine the cosmological parameters. 

Let us consider such observations from a deep survey, corresponding to supernova observations from SNAP. If we 
shun lensing, but do not know the dumpiness of the Universe, then we introduce an additional fit parameter in the 
smoothness a. While this increases the uncertainties in the other parameters, the more insidious effect is bias, or 
misestimation of the parameters. Within the context of a SUGRA model, assuming a smooth Universe while only 
utilizing lines of sight with am = 0.9 results in a bias to wq of 0.3(7, and to Wa of 0.8cr, where a is the statistical 
precision. This gives a measure of the peril of shunning lensing, and an estimate of the damage of low statistics. 

We emphasize that homogeneous observations that do not select special lines of sight do not run into such bias — it 
is better to ignore lensing than shun it, and best of all is to acknowledge and include it as done in this article. Could 
we form a probability distribution such as simulated in the Stochastic Universe Model from an analytic convolution 
of clumpy lines of sight? That is, could we incorporate both underdense and overdense lines of sight in such a way as 
to conserve flux and allow a simpler treatment than Monte Carlo analysis? 

Two difficulties exist with this method. First, overdense lines of sight reach a limit at some critical a, e.g. a > 25/24 
for the flat, pure matter model, where the increased focusing forms caustics — i.e. the light rays focus before reaching 
the observer, and the angular distance vanishes there. But this can be overcome by analytic continuation. The second 
problem is that averaging must be done not in terms of density, i.e. restoring (a) — 1, but in terms of flux conservation 
{r^'^) — rpj^yir, where FRW means a completely smooth, a = 1 Universe. This requires the use of a redshift dependent 
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smoothness parameter, which can be managed, but it also becomes a recursive relation fcf. iLinde^ l)1988(l . Appendix 
C), which is no easier to deal with than the physically motivated PDF formalism used in the text of this paper. 

Thus, our conclusions are that the Stochastic Universe Model, and the fitting functions we derive, offer the best way 
to realistically include lensing from the presence of structure in the Universe. 
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